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Abstract Owing to effects arising from quantum elec- 
trodynamics (QED), magnetohydrodynamical fast modes 
of sufficient strength will break down to form electron- 
positron pairs while traversing the magnetospheres of 
strongly magnetised neutron stars. The bulk of the en- 
ergy of the fast mode fuels the development of an electron- 
positron fireball. However, a small, but potentially ob- 
servable, fraction of the energy (~ 10 33 ergs) can gener- 
ate a non-thermal distribution of electrons and positrons 
far from the star. This paper examines the cooling and 
radiative output of these particles. Small-scale waves may 
produce only the non-thermal emission. The properties 
of this non-thermal emission in the absence of a fire- 
ball match those of the quiescent, non-thermal radia- 
tion recently observed non-thermal emission from sev- 
eral anomalous X-ray pulsars and soft-gamma repeaters. 
Initial estimates of the emission as a function of angle in- 
dicate that the non-thermal emission should be beamed 
and therefore one would expect this emission to be pulsed 
as well. According to this model the pulsation of the non- 
thermal emission should be between 90 and 180 degrees 
out of phase from the thermal emission from the stellar 
surface. 

Keywords gamma-rays: observations • pulsars: individ- 
ual SGR 1806-20, AXP 4U 0142+61, AXP IE 1841-045 • 
radiation mechanisms: non-thermal 
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1 Introduction 

Simply put magnetars are neutron stars whose magnetic 
fields dominate their emission, evolution and manifes- 
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tations. In the late 1970s and early 1980s, a fleet of 
sensitive detectors of high-energy radiation uncovered 
two new phenomena, the soft-gamma repeater and the 
anomalous x-ray pulsar. Strongly magnetized neutron 
stars provide the most compelling model for both types 
of object, and observations over the past few years indi- 
cate that these phenomena are two manifestations of the 
same type of object. Soft-gamma repeaters exhibit qui- 
escent e mission similar to that of anomalous x-ray pul- 
sars Ce.g iRothschjld et alJll994 iMurakami et al.Ml994t 
IHurlev et alJll996t Il999l). and anomalous x-ray pulsar s 
sometimes burst i|Gavriil et alJl2003 iKaspi et al.ll2003|h 
What makes magnetars a hot topic of research is the 
rich variety of physical phenomena that strong magnetic 
fields exhibit. 

This article will focus on the quiescent emission from 
these interesting objects rather than the bursts them- 
se lves (The reader may w i sh to refer to the seminal work 
of lThomn son fc Duncan! ^199&_ for details of the burst 
but may also want to look at iHevl fc HernoTdsH <|2005b) 
for an alternative). Furthermore, the article will concen- 
trate on a possible model for the recently detected non- 
thermal emission from these objects. 

In earlier work, Lars Hernquist and I analysed wave 
propagation through fields exceeding the quantum criti- 
cal value Bqed = m 2 c 3 /eh « 4.4 x 10 13 G, and de mon- 

strated circumstances under which electromagnetic (jHevl fc Hernauistl 
1998b) and some MHD waves, particularly fast modes 
l|Hevl fc Hernouistlll999T) evolve in a non-linear manner 
and eventually exhibit discontinuities similar to hydro- 
dynamic shocks, owing to vacuum polaris ation from quan- 
tum el ectrodynamics (QED). In Paper I l|Hevl fc Hernauistl 
2005b|), we developed a theory to account for bursts 
from SGRs and AXPs based on "fast-mode breakdown," 
in which the wave energy is dissipated into electron- 
positron pairs when the scale of these discontinuities be- 
comes comparable to an electron Compton wavelength. 
We showed that, under appropriate conditions, an ex- 
tended, optically thick pair-plasma fireball would result, 
radiati ng primarily in hard X-ray s and soft 7-rays. In Pa- 
per II l(Hevl fc Hernauistl l2005al) we developed the the- 



2 



Jeremy S. Heyl 



ory of the non-thermal emission from the fast-mode cas- 
cade under a series of assumptions and approximations. 
In particular, MHD fast modes of insufficient amplitude 
to generate an optically thick fireball will still dissipate 
through pair- production, seeding non-thermal emission. 
iHevl fc Hernauistl (|2005tJl estimates the Thomson opti- 
cal depth through the pair plasma created by a passing 
fast mode (their Fig. 3). If this optical depth exceeds 
unity, one would expect a fireball to form. 

In what follows, I extend our previous investigation of 
fast-mode breakdown to estimate the spectrum of non- 
thermal emission expected outside the region containing 
an optically thick fireball with particular emphasis on 
the assumptions made in Paper II and on the angular 
dependence of this emission. 



2 The Observed Non-Thermal Emission 

What would be the typical flux of t hese small-scale fast 
modes ? lThompson fc Duncan! | l996) an d lHevl fc Kulkarnil 
|)1998|) have argued that the quiescent emission of SGR 
and AXP neutron stars may be powered by the decay 
of the magnetic field. The quiescent thermal emission 
may only be a small fraction of the total energy released 
by the decay of the magnetic field. Recent observations 
of SGRs and AXPs indicate that the ther mal radiation 
may i ndeed be just the tip of the iceberg iMolkov et alJ 
120041 iMereghetti et alJl2004l iKuiper et alJl2004l) : there- 
fore, following the discussion of Paper II, the amplitude 
of the thermal and non-thermal radiation are taken to 
be comparable, and the pair cascade operates beyond 
a certain radius from the star or equivalently below a 
certain magnetic field strength (-B max )- This model has 
two parameters, B mlLX determines the position of the two 
breaks in the spectrum and the total normalisation. The 
process of fast-mode breakdown according to Paper II 
predicts a particular relationship between the location 
of the two breaks in the spectrum and particular slopes 
between and beyond the breaks. 

Figure Q shows the observed broad-band spectrum 
of several AXPs and SGRs. To compare the spectra the 
more distant objects (AXP IE 1841-045 and SGR 1806- 
20) whose hard X-ray emission was discovered with IN- 
TEGRAL have been plac ed at the distance of 4U 0142+61 
whose optical emission llHulleman et alJ l2000|) is very 
likely to be nonthermal ijOzefeOOil) . The observed spec- 
tra depart from a power law for E < 1 keV because in 
this region the thermal radiation from the neutron star 
begins to dominate the non-thermal component. 

Because the location of the two breaks in the spec- 
trum both depend on the strength of the magnetic field 
at the inner edge of the breakdown region, the presence of 
extensive nonthermal optical emission indicates that the 
non-thermal hard X-ray emission should peak at about 
30 MeV, a factor of two hundred beyond the observed 
spectrum. The best limits in this energy range are pro- 
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Fig. 1 The spectrum produced by fast-mode breakdown is 
superimposed over the observed thermal and non-thermal 
emission from several AXPs and SGRs for models that fit 
either the optical or INTEGRAL data solely and one that 
fits both sets of data. The unabsorbed optical data are from 
IHulleman et~aTl 112001) via lCteeil 120011) for AXP 4U 0142+61. 
The uppermost bla ck sym bols are the har d X-ray band are 
fromlMolkov et alJ (12001) for SGR 1806-20. ICTereghetti et all 
(120041) obtained similar results for the SGR. The middle sets 
of points in the hard X-ray data (blue i s tota l flux and 
red is pulsed flux) are from IKuiper et aTl (12001) for AXP 
IE 1841-04 5. The green squares plo t the INTEGRAL data 
rep orted bvlden Hartog et all l|200l) for AXP 4U 0142+61. 
The d en Hartog et all 112001) results are normalised using the 
observations of the Crab bv l.Tund l|19891) . We scaled the emis- 
sion from the three sources by assuming that they all lie at 
the distance of AXP 4U 0142 +61. The assumed dist ances are 
3 kpc for AXP 4U 01 42+61 llHulleman et aT]l200?l). 7.5 kpc 
for AXP IE 1841-045 llSanbonmatsu fc Helfandll9921) and 15 
kpc for SGR 1806-20 l|Molkov et al.1120041 ). 



vided by Comptel l)Kuiper et alJ 1200*6]) and appear to 
exclude the "optical model" for 4U 0142+61. 

This conclusion assumes that the sources SGR 1806- 
20 and AXP IE 1841-045 have a similar optical excess 
to 4U 0142+61. A more conservative assumption would 
be that the hard X-ray emission does not extend far 
beyond the observations from INTEGRAL with spec- 
tral breaks at about 1 MeV and 650 eV. This situation 
is somewhat natural. The fast-mode cascade is limited 
to pairs with sufficient energy to produce photons with 
E > 1 MeV that can subsequently pair produce. Lower 
energy electrons simply cool, giving the observed cool- 
ing spectrum in the hard X-rays. The total energy in 
the non-thermal emission is also reduced by a factor of 
a few. In the context of the fast-mode cascade it is dif- 
ficult to have -Ebrcak < 2mc 2 . This model is denoted 
as the "Minimal Model" in Fig. ^ because -Ebreak takes 
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on the minimal value that makes sense physically; i.e 
w 1 MeV; this is equivalent to assuming that magnetic 
field is much weaker than the quantum-critical limit in 
the pair-production region. 

The analysis reviewed here from Paper I assumed 
that the fast-mode was travelling perpendicular to the 
field and that the amplitude of the wave and the mag- 
netic field where the shock forms are independent. The 
more detailed analysis that follows shows that the loca- 
tion of shock formation depends on the product of initial 
amplitude of the wave and its wavenumber, its direction 
of propagation and the strength of the dipole compo- 
nent of the star's magnetic field. The latter two param- 
eters could potentially be determined independently of 
the emission spectrum, leaving a single free parameter to 
describe the emission, the product of initial amplitude of 
the wave and its wavenumber. 



emission are strong affected. Additionally photon split- 
ting may play a ro le. Restricting the results to £ -C 1 
yields the results of lHevl fc Hernauistl (|2005a^ . 

Specifically in the strong-field regime, neither the pri- 
mary pairs nor the secondaries can be assumed to have 
relativistic motion perpendicu lar to the magnetic field . 
The shock travels at a velocity ijHevl fe Hernauistll997bj) . 



-, «QED v 

1 Ai 

47T 



sm 



(2) 



where the function is defined in iHevl fc Hernauistl l|l997ah 
and well approximated by the fitting formula 
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from iPote khin et alJ 120041) . 

The primary pairs typically have a Lorentz factor of 



3 Spectrum 

The previous calculation of the spectrum assumed that 
the energy dissipated by the wave over a given range 
of magnetic field strengths is constant and that the dis- 
sipation only occurs below a particular magnetic field 
strength. A direct calculation the evolution of the wave 
as it travels away from the star relaxes both of these as- 
sumptions. This determines that amount of energy dumped 
into pairs at various magnetic field strengths; however, 
it does not determine how this energy is finally dissi- 
pated. It can be emitted locally and promptly as classi- 
cal synchrotron radiation or Landau transition radiation 
or elsewhere as curvature emission. 
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The typical Landau level occupied by these primary pairs 
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3.1 Weak-Field Regime 

If the magnetic field is much weaker than the quantum- 
critical limit, the spectrum of the cooling pairs at a par- 
ticular value of the magnetic field is given by 
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where E = m e c 2 £, E hrcak 



(1) 

V(40£)ande = B/S QBD 
for B <C Bqed ■ This assumes that the final generation of 
pairs emits classical synchrotron radiation as described 
in Paper II and that only the momentum perpendicular 
to the magnetic field is dissipated as synchrotron emis- 
sion. 



3.2 Strong-field Regime 

For £ ~ 0.1 the analysis of the preceding section breaks 

down. This is signaled by the fact that Eq ~ -©break, so 
Eq. JJl, for example, does not make sense. In this regime 
both the processes of pair production and synchrotron 



independent of angle, so for £ > 27 only t he ground Lan- 
dau le vel is occupied by the primary pairs. lUsov fc Melrose! 
(1995) found a limit of £ > 0.1 for pairs to be created 
of pairs in the ground Landau level. The limit found 
here is quite a bit different simply because the pair cas- 
cade discussed here is due to photons that are essen- 
tially created above the pair-production threshold in the 
shock. The direction of the photons is constant rela- 
tive to the field direction; it is their energy that in- 
creases above the pair-production threshold; they carry 
much momentum perpendicul ar to the field. On the other 
hand, lUsov fe Melrose! l)l995|) examine the pair produc- 
tion by photons travelling nearly parallel to the magnetic 
field. The momentum perpendicular to the field is small 
so the ground Landau level is preferred at much lower 
fields. 

The Landau level of the secondaries is a factor of 
sixteen smaller, so for £ > 8 only the ground level is oc- 
cupied by the secondary pairs. In these cases there is no 
prompt emission at all from the pairs. The energy of the 
initial photon goes into the rest mass of the pairs and 
their motion along the magnetic field. If n is not large, 
the emission is better described by cyclotron emission 
than by synchrotron. The emission is not cutoff by con- 
siderations of the total pair- production optical d epth as 
in the weak-field case l|Hevl fc Hernauistll2005a|) but by 
the Landau levels of t he pairs produced in the c ascade. 
In the weak-field limit, IHevl fc Hernauistl (|2005a|) found 
that the energy of the final pairs in the cascade was a 
factor of 1,700 lower than the primaries, so for £ > 0.08 
these final pairs occupy the ground Landau level and 
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the synchrotron analysis of the cascade emission must 
be modified; the limiting field for the final pairs in a 
cascade t o be in ground Landau level is similar to that 
found bv lUsov fe Melro'sel (^995) for the pulsar cascade. 
Comparing the values of Eq and inbreak yields the same 
limiting field. 



0.8 



3.3 The Total Spectrum 

The analysis of the pair cascade in the strong-field limit 
is beyond the scope of this small contribution. Instead 
I shall focus on the regime where the entire cascade is 
in the weak-field limit. This may occur for a very strong 
wave propagating near a weakly magnetised star or for 
a relatively weak wave travelling near a strongly magne- 
tised star. 

In an arbitrary magnetic field, the evolution of a fast 
mode can be parametrized by the integral of the "opacity 
to shocking" 

fcsin6* b 6 [- On i . • , 

(4-3sin 2 6>) k 2 \ (6) 
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along the path of the wave. In this expression, b is the 
amplitude of the magnetic field of the fast mode, B is 
the amplitude of the local magnetic field, 9 is the angle 
between the direction of propagation and the magnetic 
field and k is the wavenumber of the fast mode. 

Although this quantity is not strictly an opacity in 
the sense of radiative transfer, it does characterize the 
evolution of the wave in an analogous way. In particular, 
the shape of the wave depends on the optical depth for 
shocking 



Fig. 2 The opacity as a function of angle at various mag- 
netic field strengths. The curves from bottom to top are 
for £ = 0.1, 1, 10 and 100. For weak fields, the opacity 
peaks at 6 = n/2. For strong fields the opacity peaks at 

9 = cos- 1 (l/V3) w 55°. 
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where $(x) = dlnT/dx, = d^/dx and C(3) « 

1.202. Both the opacity and the energy available for 
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The shock forms at r = 1, and the power carried by the 
wave decreases as r increases above unity. 

The quantity in brackets is the gener alization of the 
functio n £F(£) where F(£) was defined bv lHevl fc Hernauistfc rompt synchrotron emission drop as the angle of prop- 
agation departs from perpendicular. 

The total optical depth is simply the integral of the 
opacity over the path of the wave. If the wave is spherical, 
the value of 6 decreases as 1/r as the wave travels away 
from the star. For simplicity the angle of the wave prop- 
agation with respect to the field is held constant; this is 
appropriate for waves travelling radially and holds ap- 
proximately for waves that have travelled several stellar 
radii from the surface. The d issipation of the energy o f 
the wave is gi ven by Eq. 72 inlHeyl Hernaui st ( 1998b) 
and Fig. 2 in lHevl fc Hernauistl <|2005bl h 

Summing over all of the magnetic fields where the 
dissipation occurs yields 

dE± , n .dP 

= -(l-cos0)— (16) 
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Fig. and Fig. |2| depict the opacity as a function of 
angle and magnetic field strength. For weak fields the 
opacity is proportional to sin 9; while in strong fields, it 
is proportional to sin0 (2 — sin 2 9) 



The definition of X (x) from lHevl fc Hernauistl ljl997a|) 
yields 
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where P is the power remaining in the wave and 9 is 
the angle that the direction of propagation of the wave 
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Fig. 3 The opacity as a function of the strength of the 
field at 8 = -n/2 in units of k with b = B. The reciprocal 
of the ordinate is the number of wavelengths that a wave 
of amplitude equal to the background field will travel before 
forming a shock. 
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Fig. 4 Total spectrum for strong-field neutron star (-Bns — 
30Bqed) as a function of angle. The initial wave has A ~ 
100 m and b ~ 2Bqed- The curves from shortest to tallest 
are for 9 = 60°, 70° ,80° and 90°. 



makes with the magnetic field. The angular dependence 
here comes from the assumption that only the compo- 
nent of the momentum perpendicular to the field yields 
prompt and local emission. Furthermore, this equation 
also assumes that the emitting pairs are produced in a 
weak-field region; they occupy high Landau levels, so the 
classical treatment of synchrotron emission holds. 

For large distances from the neutron star, dP/dB 
is a constant and the fields are weak; this in combina- 
tion with the resu l ts of § 13.11 yields the spectra given in 
iHevl fc Herncniistl |2Q05a) and shown in Fig. 0] where the 
opacity depends on angle as in Eq. © . The vertical axis 
gives the fraction of the energy of the initial wave that 
is dissipated over a factor of e « 2.7 in energy. For the 
weak-field treatment to be accurate, both the frequency 
and the amplitude of the wave must be sufficiently small 
to delay the formation of the shock until the weak-field 
regime. This results in a wave that only dissipates a small 
fraction of its energy into pairs. 

In this case the total optical depth to shocking is 
only slightly greater than unity (tqq — 1.08) so the spec- 
trum depends strongly on the angle of propagation as de- 
picted in Fig.Q] The radiation is emitted only if r > 1, 
so the radiation is emitted only for sin#rgo > 1 if the 
bulk of the opacity was in weak magnetic fields, and 
sin 6(2 — sin 2 9)rgo > 1 if it was in weak magnetic fields. 
Because the star in this case has £ = 30 but the shock 
does not form until £ < 0.1, one would expect the cutoff 
to lie between these extremes. These two estimates yield 
cutoff angles of 67° and 33°. The detailed calculations 
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Fig. 5 Total spectrum for weak-field neutron star (-Bns — 
O.IBqed) as a function of angle. The initial wave has A ~ 20/i 
and b ~ O.OIBqed- The curves from shortest to tallest are 
for 6 = 22.5° , 45° , 67.5° and 90° . 



find that the emission cuts off within about 59° of the 
field direction. 

If the field at the surface of the star is less than about 
one-tenth of -Bqed, the discussion of § 13.11 apply, re- 
gardless of where the shock forms, yielding the results 
of Fig. 
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In this case the shock may form relatively close to 
the star and a large fraction of the initial energy in 
the wave is dissipated over a narrow range of magnetic 
field strengths where the shock forms. The bulk of the 
energy may be released in the regime where dP/dB is 
not constant so the chara cteristic spectrum discussed in 
iHevl fc Hernauistl l)2005a|) does not emerge. For large an- 
gles with respect to the magnetic field, the shock forms 
further from the star and the bulk of the energy dissi- 
pates in the constant-dP/dB regime. Again the emission 
cuts off beyond a critical angle. Here T90 = 6.82522 yield- 
ing a cutoff angle of about 8° . Because the optical depth 
in the perpendicular direction is so large, the bulk of the 
angular dependence in the emission comes not from the 
opacity but from the (1 — cos0) factor in Eg. [TBI 



4 Comparison with Observations 

The model outlined in the preceding sections yields at 
least one important prediction. The hard x-ray emission 
if it is produced by fast-mode breakdown will be largest 
where the line of sight to the pulsar is perpendicular to 
the local magnetic field. This happens when the equato- 
rial regions of the pulsar face the observer. On the other 
hand, the thermal emission from the pulsar is expected 
to be strongest when the polar regions face the observer 
(e.g. IHevl fc Hernauistl Il998a|) : consequently the hard 
and soft x-ray emission should be between 90 and 180 
degrees out of phase relative to each other. 

If the magnetic field of the neutron star consists of a 
dipole whose dipole moment makes an angle (3 with the 
spin axis. Furthermore, the line of sight makes an angle 
a with the spin axis. If a + [3 > 90°, the hard emission 
peaks at a phase 

cf> = arccos (— cot a cot (3) . (17) 

relative to the main pulse of soft emission (for a+[3 > 90° 
there can be a second pulse of soft emission from the 
second polar region). Because a and (3 lie between zero 
and ninety degrees, their cotangents are positive, so <j) 
must lie between ninety and one hundred eighty degrees. 
On the other hand if a + (3 < 90°, the hard emission 
pe aks 180° out of phas e from the soft emission. 

iKuiper et alJ l|2006|) observed the x-ray emission from 
several anomalous x-ray pulsars from 0.5 to 300 keV. 
The AXP 1RXS J1708-4009 exhibits variability simi- 
lar to that described here. In the 1.3—3.9 keV-band the 
emission peaks around phase 0.25; while above 8 keV 
the emission peaks around phase 0.55 (90 degrees away). 
The peak remains at phase 0.55 to the highest energies 
measured about 300 keV. 

The AXP 4U 0142+61 shows a simila r effect albeit 
less dr amatically. At the lower energies that lKuiper et alJ 
(2006) examined (0.5—1.7 keV), the pulse profile is dou- 
ble peaked with two peaks about 160° out of phase. 
As the energy of the radiation increases to 50 keV, the 



second peak increases in amplitude and the first peak 
practically vanishes. The AXP IE 2259+586 shows sim- 
ilar trends but the variability is only detectable up to 
about 25 keV. On the other hand the pulse profile of 
AXP IE 1841-045 does not change much from 2.1 keV 
to 100 keV. This might mean that this object is an ex- 
ception or that the soft thermal emission only begins to 
dominate at lower energies. 



5 Discussion 

This paper has examined the angular dependence of the 
non-thermal emission from AXPs and SGRs in the weak- 
field regime. It is a natural extension of the model for 
bursts and non-thermal emission presented in Papers I 
and II. When a dislocation of the surface of the neu- 
tron star is sufficiently large, the resulting fast modes 
will produce sufficient pairs to make the inner magneto- 
sphere of the neutron star opaque to x-rays, generating 
a fireball. Some small but observable fraction of the en- 
ergy initially in the fast modes is dissipated outside the 
opaque region yielding a characteristic fast-mode break- 
down spectrum. A second possibility is that the crust of 
the neutron star is constantly shifting over small scales, 
generating fast modes whose breakdown is insufficient 
to produce a fireball. In this case we would associate the 
non-thermal radiation with the quiescent thermal radia- 
tion from the surface of the star. Both are produced by 
the quasi-continuous decay of the star's magnetic field. 

Observations of the AXPs and SGRs continue to sur- 
prise, as do theoretical investigations of ultramagnetised 
neutron stars. This and the preceding papers have pre- 
sented a unified model for the thermal burst emission 
and non-thermal emission from ultramagnetised neutron 
stars. The model has few underlying assumptions: mag- 
netars produce fast modes sufficient to power the non- 
thermal emission and, more rarely, the bursts, the mag- 
netic field far from from the star is approximately dipo- 
lar and quantum electrodynamics can account for the 
dynamics of pairs and photons in strong magnetic fields. 
The model for the non-thermal emission from a particu- 
lar fast mode depends only on the product of the strength 
of the wave and its wavenumber, the angle the wave prop- 
agates relative to the magnetic field direction and the 
strength of the dipole component of the field. The inten- 
sity of the emission depends strongly on the direction of 
the initial wave. Waves that travel sufficiently close to the 
direction of magnetic field do not produce non-thermal 
emission by this mechanism at all. The emission is there- 
fore beamed perpendicular to the magnetic field lines, so 
the non-thermal emission should be between ninety and 
one hundred eighty degrees (depending on the viewing 
geometry) out of phase relative to the thermal emission 
which is expected to be strongest along the magnetic 
field lines. 
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Further calculations of the fate of the pairs produced 
in the strong field regions (£ ~ 0.1) is needed to un- 
derstand the spectra produced by fast-mode breakdown 
in general; however, the beaming of the radiation does 
appear robust. Further observations can easily verify or 
falsify this model and potentially provide direct evidence 
for the ultramagnetised neutron stars that power AXPs 
and SGRs and the macroscopic manifestations of QED 
processes that account for their unique attributes. 
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